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Q\ [ ABSTRACT 

Context. X-ray spectroscopy of active galactic nuclei (AGN) offers the opportunity to directly probe the inner regions of the accretion 

1 I ' disk. Reflection of the primary continuum on the circumnuclear accreting matter produces features in the X-ray spectrum that help to 
pL^ , explore the physics and the geometiy of the innermost region, close to the central black hole. 

^J^ ■ Aims. We present the results of our analysis of average AGN XMM-Newton X-ray spectra in the Chandra Deep Field South observa- 

; ' tion (hereafter, XMM CDFS), in order to explore the Fe line features in distant AGN up to z~3.5. 

,S^ ' Methods. We computed the average X-ray spectrum of a sample of 54 AGN with spectroscopic redshifts and signal-to-noise ratio in 

CIh, the 2-12 keV rest-frame band larger than 15 in at least one EPIC camera (for a total of 100 X-ray spectra and 181623 net counts in the 

2-12 keV rest-frame band). We have taken the effects of combining spectra from sources at different redshifts and from both EPlC-pn 

and EPIC-MOS cameras into account, as well as their spectral resolution; we checked our results using thorough simulations. We 

explored the iron line components of distant AGN focusing on the narrow core which arises from material far from the central BH 

^ , and on the putative relativistic component emitted in the accretion disk. 

I I . Results. The average spectrum shows a highly significant iron feature. Estimating its equivalent width (EW) with a model-independent 

method suggests a higher EW in a broader range. The line, modelled as an unresolved Gaussian, is significant at 6.8cr and has an 

^—{ ' EW=95 ± 22 eV for the full sample. We find that our current data can be fitted equally well adding a relativistic profile to the narrow 

^ , component (in the full sample, EW=140 ± 120 eV and 67±28 eV respectively for the relativistic and narrow lines). 

(^^ , Conclusions. Thanks to the high quality of the XMM CDFS spectra and to the detailed modelling of the continuum and instrumental 

0\ ' effects, we have shown that the most distant AGN exhibit a highly significant iron emission feature. It can be modelled both with nar- 

w^ ' row and broad lines which suggest that the EW becomes higher when a broader energy range around the line centroid is considered, 

r^^ ^ provides tantalising evidence for reflection by material both very close and far away from the central engine. The EW of both features 

• . are similar to those observed in individual nearby AGN, hence they must be a widespread characteristic of AGN, since otherwise the 

nI ■ average values would be smaller than observed. 

o ■ 

fT^ ' Key words. Galaxies: Active -X-rays: galaxies 

> ; 

•'H ' 1. Introduction approach the speed of light and, furthermore, the strong gravi- 

^% tational field of the BH itself will cause a gravitational redshift 

i^ X-ray spectroscopy of active galactic nuclei (AGN) provides one eff^ect. Since t he theoretical predictions on AGN relativistic lines 

. Cp . of the most useful probes of the properties of the accretion disks ('Fabian et ^."1989*), such features have been found in the X-ray 

in the vicinity of the super massive black holes (SMBHs). In spec tra of several Seyferts, being MGC-6-30-15 the best studied 

particular, the iron K„ line is widely used to study the SMBH ^^^ ^Tanaka et al II 1 995b 

relativistic effects, therefore giving unique information about the 

location and the kinematics ofthe accreting material. Relativistic a broad Fe Ka line can be dilut ed in the continuum i f 

effects produce a line deformation when fluorescence occurs in the signal-to-noise ratio (S/N) is low (iGuainazzi et al.ll2006h . 

regions a few gravitational radii away from the central SMBH. This demonstrates the necessity of high S/N data with good 

Given that the inner radius of the accretion disk is 6Rg around a counts statistics to detect the relativistic features in the iron lines 

non-rotating black hole (BH) and 1 .23 ^ , (Rg^^) if the BH is (IG uainazzi et al.l201 iLldeTa C alle Perez et al.l2010HNandra et all 

maximally rotating dBardeen et al.ll 19721) . the velocities involved 200^. Theoretical predictions (Ballantvne (20101) supported this 

scenario, showing that the majority of relativistic iron lines have 

* Based on observations by XMM-Newton, an ESA science mission equivalent widths (EW) of just 100 eV, highlighting the impor- 

with instruments and contributions directly funded by ESA member tance of high sensitivity surveys to explore them. In cases where 

states and NASA. the individual X-ray observations have a limited S/N, statistical 
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methods like i ntegrating or averagi ng the spectra can be used for 
their analysis. iNandra et al.l (1 1997 *) used for the first time a sta- 
tistical approach to study the ASCA X-ray spectra of an AGN 
survey: they found that a broad line is present not only in the ma- 
jority of their Seyfert sample, but also in their average spectram, 
with a shape similar to those of the individual source spectra 
probing therefore that these features are common in the low red- 
shift universe. Their more recent work (Na ndra e t al. 2006 ) on 
an XMM-Newton broad line survey confirmed this result, find- 
ing that in 1/3 of their sample the iron lines are well described by 
the sum of a naiTow line likely coming from reflection in distant 
cold matter in the inner edge of the torus, and a relativistic line 
coming from the accretion disk, much closer to the SMBH. 

The deep observation with XMM-Newton of the Lockman 
Hole was studied with a statistical approach, showing a broad 
stacke d iron line, with a relativistic shape (Streblvans ka et al.l 
l2005h . Some furt her tantahsiri g evidence for a br oad iron liri e 
was also found in iBrusa et al.l d2005.) and .Civano et al.l (l2005h : 
they computed the stacked X-ray spectra of the Chandra Deep 
Fields, cautioning that their large EW and broad profile signifi- 
cantly depend on the modelling of the underlying continuum. In 
fact, the EW of the lines found in these early works were much 
higher (of the order of many 100 eV) than those found in local 
individual AGN, which casted some doubts on their reality. 

Unlike these first spectroscopic studies of the AGN surveys 
with stacking, many recent papers have successfully detected 
narrow iron lines, while the broad lines are not clearly seen. 
Examples of such detections in clude the XM M-Newton X-ray 
observations of the AXIS ( Carrer a et aDl2007h and XWAS sam- 
ples studied bv lCorral et al.i (2008) and the subsample with 0.2- 
1 2 ke V counts greater than 1 000 of the 2XMM X-ray source cat- 
alogu e (IWatson et al.l l2()08) studied bv lChaudharv et al.1 (1201 OL 
I2OI2I) . Therefore, the detection of broad Fe lines in stacked spec- 
tra of distant AGN remains controv e rsial. 

More recently, llwasawa et al.l (1201 2l) found evidence for 
ionised iron lines in the stacked spectrum of the COSMOS sam- 
ple. Broad symmetrical lines were found in the l ow luminosity 
- low redshift AGN in the deep Chandra Fields (^F alocco et alJ 
I2OI2I) . Taking the COSMOS results into account, such symmet- 
rically broadened lines were interpreted as the combination of 
a relativistic profile and several narrow lines centred at 6.4 keV 
and higher energies, from iron at a variety of ionisation states. 
Consequently the total iron line emission will appear symmetri- 
cally smoothed ( Falocco et al. 2012) . 

In this paper, we present the result of an X-ray spectral analy- 
sis of the deepest XMM-Newton survey made so far; the 3 Ms 
obse rvation centred in the Chandra Deep Field South f Comastri 
et al. l201 1'). This sui'vev enables a much more detailed spectro- 
scopic study than earlier data on this or other deep surveys and 
allows us to study the accretion at the epoch of the peak of the 
activity in the history of galaxy evolution at z~l-2 dUeda et alj 
[2003; Ebrero et al. 2009). It covers three orders of magnitude 
in luminosity (lO"*''^ - 10'*'*'^erg s"'), which allows us to study 
the behaviour of the spectral features coming from the accretion 
disk in X-rays, and the trend found between the Fe line EW with 
the luminosity ("Iwasa wa & Taniguchilll993h . 

The paper is organised as follows. The properties of the sam- 
ple are described in Sect. |2] the method used to average the 
spectra in Sect. |3j our results are illustrated in Sect. 5] and in 
Sect. [5] we discuss the implications, summarising our conclu- 
sions in Sect.|6l Throughout this paper, we adopt the cosmologi- 
cal parameters; H o=70 km s"' Mpc"', Q.,„aiter=0.3 and Qa=0.7; 
all the counts refer to the net number of counts between 2 and 
12 keV rest-frame, unless explicitly stated otherwise. We used 



xspec V . 12.5 (lArnaudll 19961) for the X-ray spectral analysis; 
the error estimates correspond to 90% confidence level for one 
interesting parameter. 

2. X-ray sample 

We present in this section the XMM-Newton observation in the 
CDFS and its spectral analysis in 12.11 the definition of the 54 
sources sample in 12.21 and the definition of the various subsam- 
ple s in 12. 31 



2.1. The X-ray data 

The XMM CDFS sui-vey is an observation of the Chandra Deep 
Field South with XMM-Newton for a total of 3.3 Ms of full ex- 
posure time. Exposure times from the EPIC pn and MOS cam- 
eras are 2.5 Ms and 2.8 Ms respectively. The bulk of the ob- 
servations were performed between July 2008 and March 2010, 
and have been combined with archival data taken between July 
2001 and January 2002. The details on data reduction and anal- 
ysis and the source detection algorithm are presented in Ranalli 
et al. (.201 3h . 

In this work, we used the high quality X-ray spectra for 
sources detected with a significance above 8 sigma in the 2-10 
keV source catalogue, corresponding to a flux limit of ~ 2 x 
10^'^erg s^'cm"^. 

We looked for spectroscopic identifications of these sources 
in the hterature making use of several spectroscopic campaigns 



(Le Fevre et al. 2004; Szokolv et al. 2004; Le Fevre et al. 2005 
Mignoli et al 2005; Ravikumar et al. 2007.;.Kriek et al...2008l ; 
Vanzella et allboOsl-lTavlor et al.ll2009l;lTreister et al.ll20q9L van 



der Wei et a lj2005?Balestra etal."2010|; [Silverman et al.l l2010l; 
Cas ev et al] |2011; Cooper et al. 2011]). To avoid possible spu- 
rious broad line detections caused by errors in the redshift esti- 
mate, which could produce a false line broadening, we used only 
sources with spectroscopic redshifts defined as secure in the lit- 
erature. 

To characterise the spectra and to estimate their intrinsic lu- 
minosities, after grouping the XMM-Newton EPIC MOS and 
PN spectra at 20 counts per bin, the spectra were fitted with 
xspec using the;^^^ statistics between 1 and 12 keV rest-frame in 
order to fit the broadest possible band avoiding the background- 
dominated E>12 keV rest-frame band (for this spectral region 
our grouping of 20 counts per bin would not probably be a suit- 
able choice due to the background); we do not consider the soft 
X-ray band (E<lkeV) due to the onset of complex spectral fea- 
tures (warm absorption and soft excess). We used a simple pow- 
erlaw modified by Galac tic absorption, fro zen at a column den- 
sity ~ 8 X lO'** cm"2 (IStark et al.lll992h . and variable intrin- 
sic absorption. We analysed the background spectrum and found 
an emission feature coincident with the position of the Al line 
characterising many background spectra at observed-frame en- 
ergy ~ 1.5 keV (.Carter & Read.2007.) . However, this feature is 
not expected to influence significantly our iron line analysis, as 
there is only one object at z higher than 3.27, value which would 
shift the Al feature to 6.4 keV. To assess whether this system- 
atic feature influences the output of our fits, we repeated the fits 
removing the channels corresponding to 1.5 + 0.1 keV. The fit 
results show that the presence of this background feature does 
not really influence the output fit parameters very much; com- 
paring for example the intrinsic absorption estimated in X-rays, 
we found that the difference is negligible in the majority of the 
sample, with the exception of a handful of spectra (seven, for 
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Table 1. Properties of the full sample with S/N> 15 and its subsamples 



Sample 


A'.,, 


N,, 


^-12 


A'5-8 


<z> 


<log(L)> 


<A'h,22> 


A'h.22 


{S/N) 


a 


£6 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(ergs-') 
(7) 


(10-^ cm--) 
(8) 


(10^- cm-^) 
(9) 


(10) 


(11) 


(eV) 
(12) 


Full 


100 


54 


181623 


40863 


1.34 


43.74 


1.48 


0.08 


29.8 


0.35 


119 


Unabs 
Abs 


52 
29 


32 
17 


110523 
45896 


21071 
14373 


1.27 
1.28 


43.78 
43.64 


0.02 
4.90 


0.0 
2.41 


33.60 
24.88 


0.34 
0.37 


122 
116 


z< 1 
z> 1 


39 
61 


20 

34 


72672 
108951 


18148 

22715 


0.64 
1.80 


43.21 
44.07 


1.61 
1.39 


0.17 
0.04 


30.0 
28.8 


0.35 
0.35 


95 

134 


Sey 
QSO 


61 
39 


33 
21 


91511 
90111 


20987 
19876 


0.91 
2.02 


43.35 
44.34 


1.35 
1.68 


0.11 
0.04 


26.2 
34.1 


0.34 
0.36 


108 
138 



Notes. Unabs and Abs subsamples have and log(A'H) < 21.5 and log(A'H) > 21.5 respectively; Sey and QSO subsample have log(L) < 44 and 
log (L) > 44 respectively (see text). Columns: (1) (Sub)Sample; (2) Number of spectra; (3) Number of sources; (4) number of net counts in 2-12 
keV rest-frame; (5) number of net counts in 5-8 keV rest-frame; (6) average redshift; (7) Logarithm of the average rest-frame 2-10 keV luminosity 
in units of erg s-', corrected for Galactic and intrinsic absorption; (8) and (9) average and median of the intrinsic column density in 10^^ cm-^; (10) 
average S/N of the (sub)sample between 2 and 12 keV rest-frame; (11) slope of the power-law dependence between the width of the simulated 
unresolved lines and their central energies E = E^" (see Sect. 3.3); (12) width at 6 keV of the simulated iron lines (see Sect. 3.3) 
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Fig. 1. Spectral signal-to-noise-ratio between 2-12 keV of the Fig. 2. Redshift distribution of the full (54 X-ray sources) sam- 
fuU (100 spectra) sample. pie. 



which we used the fit parameters obtained ignoring these chan- 
nels). The luminosities were estimated from the best fit models 
of each spectrum, corrected for Galactic and intrinsic absorp- 
tion, and calculated between 2 and 10 keV rest-frame. These 
are the luminosities used in the remainder of this paper for each 
spectrum. In the case of sources with both EPIC MOS and pn 
spectra, when a single value per source is presented (e.g. Fig.|3]l, 
we have used the value corresponding to the spectrum with the 
highest S/N. 

In the remainder of the paper, the luminosities are always 
indicated in erg s-' and the intrinsic absorption will always be in 
units of cm--. 



2.2. Sample definition and properties 

It has been found empirically that a good S/N is essential for a 
significant detection of a broad line with the EPIC cameras (de 
la Calle Perez et al. l2010l:lGuainazzi et al.ll201 ll) . 

To maximise the average S/N of the sample, taking the level 
of the background of the XMM-Newton observation into ac- 
count, we limited our investigation to the 100 spectra coiTe- 
sponding to 54 unique sources (49 spectra from the EPIC-MOS 
cameras and 51 spectra from the EPIC-pn camera) characterised 
with S/N>15 estimated between 2 and 12 keV rest-frame (the 



band used in the analysis of the average spectra described be- 
low). This selection of the sample represents a compromise be- 
tween the need to avoid including noisy sources and the neces- 
sity to have a high enough number of sources in the sample for 
good statistics. We discarded the two brightest sources whose 
number of counts (about 20000 net counts each) is significantly 
higher than those characterising the other sources (each of the 
sources analysed here have less than 10000 net counts in the 2- 
12 keV rest-frame band and the total accumulated counts by the 
54 sources is about 180000). In this way we avoid the estimated 
average spectra being dominated by the two brightest sources, 
which would be contributing 10% of the total counts e ach. These 
sources are: PID 319 and PID 203 ('R ananietall2013h . the latter 
shows a remarkable and variable iron line. 

The high quality of the spectra can be seen in the Fig. [1] 
where the distribution of the S/N (in 2-12 keV) of the 100 se- 
lected spectra is shown. We can see in Fig.|2]the distribution of 
the spectroscopic redshifts of the sources included in the sam- 
ple: 21 AGN in our sample are located at redshifts z~l-2 and a 
dozen AGN at earlier epochs {z ~ 2 - 3). Note, from this dis- 
tribution, that the rest-frame 2-12 keV band used in this paper, 
coiTesponds to an observed band of 1-6 keV at a typical z~l, 
therefore avoiding strong fluorescent signatures in the detectors' 
background. 
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Fig. 3. Top panel: Distribution in the Luminosity -redshift plane 
of the 54 XMM CDFS sources (pink circles) and of the 33 
Chand ra Deep Field South sources (blue stars) by ' Falocco et al.l 
( I2OI2I) . Bottom panel: Distribution in the Luminosity-redshift 
plane of the 17 absorbed (cyan stars) and 32 unabsorbed (pink 
circles) subsamples. 



Figure [3] shows the distribution of the sample sources in the 
luminosity - redshift plane. In this figure, only one entry per 
physical source is presented: in sources where both the spec- 
tra from MOS and from pn are available, we display the result 
for the highest S/N spectrum. The CDFS observation by XMM- 
Newton allows a better coverage of the luminosity-redshift plane 
than our previous work on the Chandra spectra of the sources de- 
tected in the deep Chandra fields (see top panel of Fig.[3]l. The 
XMM CDFS survey is able to characterise with the highest sig- 
nificance reached so far the properties of distant X-ray selected 
AGN over three orders of magnitude of continuum luminosity 
and a broad span of redshifts, up to ~ 3.6. Despite the good 
counting statistics of the individual spectra, the kind of analysis 
of the iron line we carry out in this paper requires much higher 
statistics, hence we perform stacking. 

2.3. Subsample definition 

The sample characteristics and properties are listed in Table 1. 
We constructed two luminosity subsamples with similar number 
of counts in each, resulting in a Seyfert subsample with log (L) < 
44 (hereafter, we call refer to this subsample as Sey) and a quasi 
stellar object (QSO) subsample with log (L) > 44 (hereafter, we 
call this subsample QSO). 



Fig. 4. log (A^h) distribution of the full sample. A^h is in units of 
10^^ cm"^; all the values log (A^h) < 20 have been replaced with 
log(A^H) - 20 for plotting purposes only. 



We also segregated the full sample in two redshift subsam- 
ples with the same number of counts. The low-z subsample in- 
cludes all the sources with z < 1 . The high-z subsample includes 
all the sources with z> 1. We did not define more than two sub- 
samples because we preferred to have as many source counts in 
each to attempt a sensitive assessment of the iron line trend with 
these parameters. 

Finally, we separated the full sample in two subsamples of 
different intrinsic column density estimated in X-rays (We sub- 
samples') setting the threshold at log (A^h) = 21.5. This is often 
used as the threshold below which the interstellar medium of the 
galaxies hosting the AGN could provide the obscuration of broad 
lines in optical spectra (Caccianiga et al. 2007). The A^h distri- 
bution (see Fig. nil extends from log (A^h) < 20 to log (A^h) ~ 23. 
The Nh shown in the figure have been obtained from the fits to 
the individual spectra. For sources with both MOS and pn spec- 
trum, we present the result from the best S/N spectrum. We can 
also appreciate in this Figure that around log (A^h) ~ 21 .5 there is 
a clear gap between lower and higher absorption sources, mak- 
ing this value a 'natural' choice for our sample. 

The two A^H subsamples have been constructed consider- 
ing the errors ANn on the A^h obtained in the fits to the in- 
dividual spectra. The sample, which includes the sources for 
which at 90 % of confidence log(A^H + AA^h) < 21.5 (here- 
after called 'Unabs' subsample) accumulates about 111000 net 
counts in the 2-12 keV rest-frame band, while the sample with 
log (A^H - AA^h) > 21.5 at 90 % of confidence (called hereafter 
'Abs' subsample) contains about 46000 net counts in the same 
band. Nineteen spectra have not been included in the two above 
mentioned subsamples (i.e., the 90% error bars around the best 
fit A^H values cross the log (A^h) = 21.5 borderline that separates 
'Abs' from 'Unabs' sources). These unclassified spectra contain 
25204 and 5419 counts respectively in the 2-12 keV and the 5-8 
keV rest- frame bands. 

We note that, despite its name, the 'Abs' subsample con- 
sists mostly of moderately absorbed AGN (median A^h ~ 2.4 x 
lO^^cm"^). Our sample does not contain any Compton thick 
sources, so we do not expect to detect in the average spec- 
tra of the absorbed sample any strong absorption feature in the 
iron line region. The two A^h subsamples are characterised by 
a similar distribution in the luminosity-redshift plane (see bot- 
tom panel of Fig.|3]l, allowing to assess whether the intrinsic ab- 
sorption alone can influence the detection of the X-ray spectral 
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features of the AGN. In the Fig. [3] one entry per each physi- 
cal source is shown (if the spectra from both MOS and pn are 
available, we plot the best S/N spectrum). In the Fig. [3] the 54 
XMM CDFS sources have higher luminosities at the same red- 
shift range than the 33 Chandra CDFS AGN. According to the 
Iwasawa-Taniguchi effect, a higher luminosity corresponds to a 
lower iron line EW, thus the spectra from the XMM CDFS will 
not be particularly advantaged in the iron feature detection. 

3. Analysis method 

The aver aging method used here was presented for the first time 
dCorral et al...2008..201 1.) in the study of X MM-Newton spectra 
of the XMS-XWAS surveys dCarrera et all 120071: iMateos et all 
1200 8*) and for the XBS surveyTbella Ceca et al.''2004'); later, it 
was adapted and tested for the deep Chandra Fields bv Falocco 
et al. (l2012h . In this work we refined the technique of analysis 
introducing new methodologies that will be described in Sect. 
[33]andl4l] 



3.1. Averaging method 

After fitting the spectra of every individual source as described in 
Sect. 12. II we read into xspec again each unbinned, background- 
subtracted spectrum, and we used the corresponding best-fit (to 
the binned spectrum) model to extract the unfolded spectrum 
compensating for the effective area of the instrument (euf spec 
command in xspec). The unfolded spectra depend on the model 
used for this operation. In Falocco et al. (2012), we have checked 
the effect of the correction for detector response on the average 
spectrum, demonstrating the independence on the model used 
in the unfolding process above ~2 keV, even for models clearly 
different from the input spectrum. Our simulations (as explained 
in Sects. 13.21 and I3.3l l take any effect of the uncertainty in the 
continuum shape into account, since each simulated spectrum is 
fitted and the corresponding best fit model is used for the unfold- 
ing, instead of the input shape. 

Once we obtained all the un folded spectra, i n analogy to 
what was done in Falocco et al. (2012) and in ICorral et al.l 
(l2008l) . we applied the following procedure: 

1 . Correction for Galactic absorption, 

2. Shift the spectra to rest-frame, 

3. Re-normalise each spectrum by dividing by its integrated 
flux in the 2-5, 8-10 keV spectral range 

We used the 2-10 keV rest-frame band in the last step (ex- 
cluding the 5-8 keV rest-frame range were the putative iron line 
is expected to contribute most) to avoid giving too much weight 
to absorbed sources, which are in much higher pr oportion of the 
full sa mple with respect to our previous work in JFalocco et al.l 
( 1201 2). where we used only the 2-5 keV rest-frame band for the 
normalisation. 

To minimise the effect of the Al spike at 1 .5 keV in the spec- 
tra from XMM CDFS (discussed in Sect. 12.11 1. we ignored the 
bins corresponding to rest-frame energy in the range 1 .4- 1 . 6 ke V 
observed frame (compensating for the lost flux by interpolating 
linearly to the full band) when estimating the continuum for the 
normalisation and also in the final average of both the real and 
the simulated data. 

We binned the spectra using an energy grid with 100 eV-wide 
bins and we finally averaged them. The errors were estimated as 
the dispersion of the fluxes around the average. 

We have found that if we use instead error bars obtained 
from error propagation of the individual source bins, as we did 



in JFalocco et al.l (l2012h . the combined error bars were smaller, 
resulting in very high reduced ;i'^ values. We believe that this is 
because, with these new data with more net counts, the intrinsic 
spectral dispersion dominates over the counting noise. 

3.2. Simulations of the continuum 

We have estimated the expected shape of the underlying contin- 
uum in our spectra using comprehensive simulations. 

We simulated the spectrum of each source 110 times using 
the best-fit parameters of the simple absorbed power law model 
from the individual fits to the real spectra (see Sect. 2.1). To each 
of these 110 simulated samples we applied the same method as 
the one used for the observed sample (spectral fitting, unfolding 
to correct for the instrument's energy-dependent response, cor- 
recting for Galactic absorption, shift to rest-frame, normalising, 
averaging and calculation of errors from the dispersion). After 
this, we adopted as the underlying continuum the median of the 
simulated continua. The median is preferred over a simple aver- 
age because it is a more robust estimate of the central value of a 
distribution. 

We present in Figs. |5] |6] |7] our average spectra together 
with their corresponding simulated continua. In those figures we 
present: the average observed spectrum (data points), our sim- 
ulated continuum (continuous line) and the Icr confidence lines 
of our continuum simulations (dashed lines). The latter were es- 
timated from the 68% percentiles around the median. 

An intense iron line characterises the observed spectra of the 
sample and all the subsamples that we constructed. Its profile 
appears visually broader than the instrumental resolution, in the 
full sample (shown in the top panel of Fig.|5]|. 

3.3. Simulations of the Fe Line 

The spectral resolution of the X-ray detectors and the averag- 
ing process can widen any narrow spectral features in the X- 
ray spectra making them to appear artificially broader than they 
actually are. The resulting effective spectral dispersion changes 
with energy and, consequently, our ability to resolve spectral 
features will change along the spectrum. To properly study the 
iron line profile in the XMM CDFS, we need to estimate our ef- 
fective spectral dispersion in the energy band we are analysing. 
To do this, we ran high S/N simulations of unresolved (cr - 0) 
Gaussian emission lines centred at rest-frame energies between 1 
and 10 keV, in steps of 1 keV with input EW of 200 eV (one sim- 
ulation for each source -at its observed redshift-, and for each en- 
ergy). The input continuum model was a powerlaw with F = 1.9 
without any intrinsic absorption. We corrected the spectra for 
detector response (model hng them with a p ower law with pho- 
ton i ndex 1.9) as done in lCorral et al.l (l2008h and lFalocco et al.l 
(l2012h : then we applied the same treatment described in Sect. [3] 
We then studied how the width of the simulated lines changes 
wi th the central energy. In addition, in analogy with what we did 
in JFalocco et al.l (l2012h . we simulated lines centred at 6.4 keV 
and 6.9 keV to estimate the resolution at the peak energy of the 
line from neutral and ionised iron. We recovered the input pa- 
rameters of the continuum very well, reproducing a powerlaw 
continuum with F ~1.9, while the Gaussian lines appear broad- 
ened by S ~ 100 eV. A similar value of the spectral resolution 
given by the stacking procedure and the X-ray instruments was 
found in the study of the XMM-Newt on X-ray spectra in the 
COSMOS sample dlwasawa et al.ll2012h . 
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Fig. 5. Average observed spectrum (data points) with its simu- 
lated continuum (continuous line) and the Icr confidence levels 
of our simulations (dashed lines). From top to bottom: full sam- 
ple, Unabs subsample, Abs subsample. The dotted line in the top 
panel represents the intrinsic resolution at 6.4 ke V obtained from 
the simulations of the iron line described in Sect. 13.31 



We studied the width of the average simulated lines as a 
function of their centroid energies: the result shows, for the full 



xO.34 



sample, S = 119 eV I g^y) ■ ^^^ values of the Z obtained in 
this work include the contributi ons from the spectra l dispersion 
of the instruments, E ~ 55 eV ( Haberl et al.ll2002l) . as well as 
the widening introduced by the averaging method. The parame- 
ters of this !,(£) dependence (the slope and the S at 6 keV) for 
the full sample and its subsamples are in Table |2] This energy- 
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Fig. 6. Average obsei-ved spectrum with its simulated continuum 
and the Icr confidence lines. Top panel: z<l subsample, bottom 
panel: z>l subsample. Symbols and line styles as in Fig.|5] 



dependent instrumental response has been taken into account in 
the spectral fitting of the average spectra. 

We emphasise that this limited spectral resolution essentially 
means that any features narrower than about 0. 1 ke V are very 
likely due to statistical fluctuations. 



4. Results 

We discuss in this section the results of our analysis of the av- 
erage XMM-Newton X-ray spectra of several subsamples in the 
XMM CDFS, estimating in particular the iron line significance 
and shape. A new way to estimate the line EW independent of 
the fits to the average spectra has been applied and it is described 
in 14.11 We describe the detailed spectral fits to the data in the 
Sect. 14:2] 



4.1. Model-independent estimate of the iine EW 

We compared the features found in our average obsei-ved spec- 
trum with the average simulated spectra to estimate the line 
EW (see Table |2]i. To do this, we used the following energy 
ranges: 5.5-7 keV (that we call, for convention, bulk), 6-6.8 
keV (called neutral), 6.2-6.6 keV (called narrow), 6.8-7 keV 
(called ionised), chosen to represent several spectral features 
from relativistically-broadened neutral Fe-K lines (5.5-7. keV), 
and from narrow neutral Fe lines (6.2-6.6 keV), to symmetrically 
broadened lines (6.0-6.8 keV) to lines from H-like Fe (6.8-7.0 
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Fig. 7. Average observed spectrum with its simulated continuum 
and the Icr confidence Hnes. Top panel: Sey subsample, bottom 
panel: QSO subsample. Symbols and line styles as in Fig.|5] 



4.2. Spectral fits to the full sample and subsamples 

To characterise the spectral features of our average spectra and 
to understand their physical origin, we fitted them in the full rest- 
frame 2-12 keV band (ignoring the E<2 keV rest-frame band in 
this analysis because we are not interested here in studying the 
soft excess). 

In our analysis of the deep Chandra fields of iFalocco et al.l 
(1201 2l) the average simulated spectrum was employed as a con- 
tinuum in the spectral fits. That approach does not take the un- 
certainty in the shape of the underlying continuum into account. 
In the current work, we preferred to construct an empirical con- 
tinuum with a flexible, smooth shape. To take the instrumental 
line broadening discussed in Sect. l3.3l into account, we smoothed 
the input spectral model with a Gaussian convolution model 
(gsmooth in xspec). This model smooths the spectral features 
with a Gaussian with a variable sigma, which changes as a func- 
tion of the energy (the powerlaw slope and the 2 found in our 
simulations of the unresolved lines in Sect. I3.3l and summarised 
in Table|2]are used), hence, all the parameter values given below 
are 'intrinsic' values, i.e., with the instrumental/method resolu- 
tion already subtracted. 

Our continuum is an absorbed powerlaw plus a Compton 
reflection component (both represented by a pexrav model in 
xspec), fitted individually to the average spectrum of the full 
sample and each subsample, excluding initially the 5-7.2 range 
(to avoid the iron feature region). 

We reintroduced the channels between 5 and 7.2 keV and 
we show in Col. 2 of Table |3]the corresponding ;^'^. The simple 
continuum model fits the overall continuum reasonably well in 
all cases, as the values of ;^'J indicate. We found an excess in the 
iron line region that we model with a simple phenomenological 
Gaussian (see 14. 2. Il l and with a more realistic model (see l4.2.2"T l. 
The significance of any fit improvement is calculated, from the 
X^, using P(Av,Ax^) = P(Av/2, A y^/2), where P(A v/2, Aa^^/2) 
is the incomplete gamma function (iPress et al.lll992h . 



keV), taking our spectral resolution (E ~ 0.1 keV) into account. 
We did not considered an intei-val centred at 6.7 keV to repre- 
sent the line from He-like Fe because the average spectra do not 
visually display this feature (see Figs. 6-8). 

We estimated the EW of the iron line using the 110 
simulated continua and used the same energy bands de- 
scribed above. For each simulation we calculated the integral: 

EW=j^^' dE'^^^^^ ~ 2 A(£,)^^^^^^ where A(£,) is the 
width of the bin with centre £,; TiEi) represents the average 
observed spectrum, C(Ej) the average continuum. The values of 
the median of the EW for each subsample and energy range are 
shown in Table|2] The confidence levels in that Table are calcu- 
lated from the 68% percentiles around the corresponding medi- 
ans. 

There is evidence for a significant decrease of the EW of 
the narrow line with increasing average luminosity, in agree- 
ment with the 'Iwasawa-Taniguchi effect' (llwasawa & Taniguchil 
Il993h . No significant trend with redshift has been found. The 
EW estimated in the 'bulk' range is higher for the Unabs sub- 
sample than for the Abs subsample, as expected (because unab- 
sorbed AGN are seen preferentially face-on), we should note that 
the Nh values covered by the absorbed sample that we defined 
in this work are too low to affect significantly the continuum 
above 6 keV, thus the iron line region. The same consideration 
should be applied to understand also the lower EW estimated in 
the 'neutral range' for the Abs subsample, despite a higher EW 
of the narrow core expected when absorption is significant. 



4.2.1. Iron line fitting with phenomenological models 

After having re-introduced the bins between 5 and 7.2 keV, that 
contain the iron line, we wanted to model the excess found in the 
iron line region. The first model we used is a simple phenomeno- 
logical Gaussian. 

We set the parameters of the Gaussian as follows: 

- Fixed cr-O and fixed centroid energy at 6.4 keV: this allows 
us to estimate the significance of the narrow component of a 
neutral Fe K line (hereafter 'fi/aw - o-fixe/, in the first line 
of each sample in Table|3]l 

- Fixed centroid energy at 6.4 keV and free cr: it studies the 
significance of a possible broad component and estimates its 
width (hereafter 'E fixed - cr/ree', in the second line of each 
sample in Tabled 

- Free centroid energy and fixed cr-Q: it considers the presence 
of a narrow Fe component centred at energies >6.4 keV and 
estimates its centroid energy (hereafter' E free -cr/ixe/, in the 
third line of each sample in Table O 

- Free centroid energy and free cr: considers both iron emis- 
sion centred at energies > 6.4 keV and broad line emission 
(hereafter ' E free - cr free', in the fourth line of each sample in 
Table EJ 

In Col. 4 of Table [3] the probability of any fit improvement 
with <Jfixed-E fixed with rcspcct to the fit with the continuum only 
is presented (in the first line of each sample). The probability of 
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Table 2. Median of the EW calculated in different energy ranges using the simulations. 



Energy band Bulk Neutral Narrow Ionised 
sample 5.5-7. 6.0-6.8 6.2-6.6 6.8-7. 



Full 



Unabs 






Abs 



z<l 



Z>1 



Sey 



252%, 

254+^;'; 

176+^''' 

-— 4?- 

283% 



160%^ 



194- 



117+1^ 
"156^ 



150+ 



166% 



865r 



160%° 

213%? 



118% 



148+ 



QSO 



71 



TW" 



87+ 



^ 



"66!l3 



-^ 



26+ 



38%^ 
1 c+K 



37+ 



21 + 16 

_ij-i7_ 
47+io 



<32 



Notes. The errors include both statistical errors propagated through the stacking process as well as the uncertainties in the continuum from the 1 10 
simulations. The range is expressed in keV and the EW in eV. 



the improvement of Efree - (T fixed, Efree - CTflxed and Efr^e - 

CTfree with respect to cr fixed - E fixed model is shown instead in 
the second, third, fourth line corresponding to each sample. 

The cr given in Table |3] is intrinsic, since our effective spec- 
tral resolution has already been taken into account through the 
convolution model g smooth. 

The values of the photon index of the average spectra are 
rather flat in all cases with the exception of the Unabs subsam- 
ple. Fitting the simulated continuum of the full sample with the 
model gsmoothCpha '' pex) we found continuum parameters 
consistent with those reported in Table|3](r - 1 .66 + 0.01 for the 
simulated continuum). 

The values of the column density Nh obtained in the fits to 
the average observed spectra are moderate (see Table |3]and Fig. 
4 for details), recovering the moderate absorption in individual 
spectra (see Table |2]l. The significance of the emission line at 
6.4 keV is 6.8cr in the full sample and always higher than 3.9 
cr in its subsamples (Col. 4 of Table [3]l. Leaving the Gaussian 
width and/or the central energy free in the spectral fitting is not 
required by our data in our full sample and the majority of the 
subsamples that we constructed. Given that the fit with free cen- 
troid energies gives output central energies equal to 6.4 keV, we 
do not add in the following analysis a further line from ionised 
iron. 

The QSO subsample has a worse statistical description with 
respect to the other subsamples (see the x^/dof in Table |3]l, 
due to relatively strong residuals above about 8-9 keV rest- 
frame which cannot be accounted for by the actual modelling. 
Moreover, the QSO subsample is characterised by a relatively 
low spectral dispersion and thus flux errors significantly smaller 
than the other subsamples. 

4.2.2. Fits with pliysical models 

Having established in the previous Sections the presence of Fe 
K emission features in the average spectra of our full sample 
and subsamples, we will now use physically motivated models 
to characterise its shape. The iron line model is overlaid to the 
same continuum model considered in the previous section; its 
profile is modelled by a Gaussian with fixed energy at 6.4 keV 
and no intrinsic width described in the previous section (that rep- 
resents fluorescence from material far from the central BH) and 
a laor component (representing fluorescence from the accre- 
tion disk arou nd a rotati ng BH). The laor model used here was 
developed bv iLaorl (11991) to test the accretion disk hypothesis 
and it assumes a Kerr metric, valid in an accretion disk around 
a rotating BH. Several other models are available in the liter- 
ature but we chose to fit our average iron lines with one of the 



simplest ones because, as shown below, our data clearly does not 
warrant a more sophisticated modelling. We set the laor param- 
eters making simple considerations about the geometry and the 
physics expected in our sources. 

The model adopted assumes emission from an accretion disk 
with emissivity index -3, in the hypothesis that the continuum 
source is located in the disk axis and at low height from the disk 
surface. Since our sample is dominated by unabsorbed sources 
that constitute 59 % of the full sample (see Table |2]i, we expect 
that the accretion disks are seen preferentially face-on, under the 
hypothesis that the accretion disk and the torus axes are aligned. 
The accretion disk is assumed to extend from 1 .23 to 400 R^ as it 
is the case for a disk around a maximally rotating BH. Hereafter, 
we will refer to the continuum plus the fluorescent lines (both 
narrow Gaussian and laor) as to the 'two-component model'. 

The fits to our data and the ratios between data and model 
are shown in Figs. from|8]to[T0] the corresponding fit results are 
given in Table |4] In the figures, a dotted vertical line has been 
placed at 6.4 keV energy for reference. In Fig. |8] we present, 
for comparison, the fits to the narrow Gaussian centred at 6.4 
keV (the E fixed - cr fixed described in the previous paragraph) in 
the top panel and the fits to the two-component model in the 
second panel. The fit is visually very similar in both cases. This 
impression is confirmed by looking at the first row of Table^that 
shows that, although the fit is better including both components, 
the improvement is statistically significant only at 2.1cr. 

The formal 90% confidence levels on the EW of both compo- 
nents almost reach zero in all cases, this is a consequence of the 
very strong coupling between their fluxes, as will be discussed 
below (see also Fig. 112b. 

We show in the last two columns of Table|4]the;t'^/v and the EW 
of the relativistic component obtained after setting the flux of the 
narrow component to zero. 



4.2.3. Fits witli compreliensive models 

We finally made a test using a full-inclusive model that considers 
the entire set of spectral components expected in AGNs: the iron 
K„ at 6.4 keV, the Fe K/^ at 7.05 keV, the Ni K^ at 7.47 keV, the 
Compton shoulder to the Fe K^, (approximated to be a Gaussian 
at 6.3 keV with width 35 eV) and the Compton reflection compo- 
nent of the continuum. This last component describes reflection 
in a slab geometr y. The full-inc l usive model, namely pexmon, 
was presented by iNandra et al.l (l2007h . We used two pexmon 
components to consider reprocessing from material close to the 
central engine and far away from it. If the emission comes from 
the innermost regions, the pexmon will be smeared by the rel- 
ativistic effects due to the proximity of the central SMBH. To 



S. Falocco et al.: Iron lines in the XMM CDFS deep survey 



data and folded model 



Full sample 




Energy (keV) 
data and folded model 







Full sample 








> 


0.05 


^-^ 


^^ 




1 

i 

■D 

to 
E 


0.02 

0.01 

5x10-= 
1.4 






"^'"%^ 




- 




|- 


5 


1.2 

1 

0.8 


I+++I 1 , i + i ' ■ 


...++4., 


f,++. ,i+i„w4l*iM„lil,iH,'ll 1 ,1 


1 1 




1 ^T-f+j+Tf+^.+T 


■■VTyTni-tt|T'.ti'p|i 


-ill 



Energy (keV) 



data and folded model 



Unabs subsample 




1.2 

1 
0.8 



.++ 1 ^ .+++V.t i l i tt iii^ii.iiiiittin - - 



1 



Energy (keV) 
data and folded model 



Abs subsample 



T3 0.01 - 

E 

° 5x10-^ 

1.4 F" 
1.2 

0.8 
0.6 




l|l+i i + i HA ' A#y 




Energy (keV) 



Fig. 8. Full sample. Top panel: Fit with Efjxed - cr fixed- Bottom Fig. 9. Fit with the 'two-component model'. Top panel: fit of the 
panel: fit with the 'two-component model' . Unabs subsample. Bottom panel: fit to the Abs subsample. 



account for these effects we convolved one pexmon component 
with the relativistic kernel kdblur2 which smears both the flu- 
orescent lines and the Compton reflection. The resulting model 
will be absorbed by circumnuclear material (pha model). The 
second pexmon component, instead, will not be smeared or ab- 
sorbed because it originates in external regions. 

The conv olution model kdbl ur2 used here has a similar set- 
up to that of Nandra et al. (2007) to allow a useful comparison 
with iron lines detected in individual spectra of AGN with high 
spectral quality. The model considers two emissivity laws, with 
slopes of and 3, and a break radius which defines the border be- 
tween the two regions initially fixed to 20 Rg for simplicity. We 
considered an accretion disk around a maximally rotating BH, 
extended from 1.235 to 400 Rg, to allow a comparison with the 
results presented in the previous paragraph. We then thawed the 
inner radius and the break radius, but the results do not change 
significantly: in the first case we found an inner radius of < 62Rg 
witiix^ /dof - 83.9/93 ; in the second case we obtained a break 
radius of < UlRg With x^ldof = 84.1/93. 

The comprehensive model just described allows a descrip- 
tion of the continuum more accurate than the phenomenological 
and the 'two component' model. The values of the photon in- 
dex increase from 1.63;^^^^ reported in Table|3]to the 1.71+0.06 
value of Table [Sj although they remain compatible within their 
error bars. The reflection factors and their 90% confidence lim- 
its, -0.35±0.18 and 0.45^035, are both different from zero, sug- 
gesting that the two components have been found (see Sect. l5.4l i. 



5. Discussion 

We discuss in this Section the results obtained with the simula- 
tions, with the fits to a phenomenological model, and with the 
fits to physicaUy-motivated models. 

5.1. Discussion based on tlie simulations 

Accurate modelling of the continuum is critical in determining 
the iron line properties of AGN and for this reason we devel- 
oped a new method to do it independently of the spectral fitting. 
The EW discussed here are estimated with respect to a contin- 
uum integrated between 5.5 and 7. keV in the full sample (see 
Sect. 14. 11 1. We expect at least two Fe line components within 5.5 
and 7.0 keV: a narrow Fe line from neutral material with EW 
of 50-100 eV ubiqui t ously se en both in nearby and distant AGN 
dNandra et al.ll2007l: ICorral et al...2008. .2011.) and a broad line 
component with a distribution spanning from virtually zero to 
3 00 eV an d a median around 170 eV, as in the Seyferts sample 
of ldelaC alle Perez et al. (2010). The EW of 203+i^j^ eV found in 
the full sample is therefore entirely consistent with having both 
components and, most importantly, not consistent with having 
only a narrow component (which is invariably smaller than 100 
eV in AGN at low redshift). The estimates of the total EW in the 
subsamples also are compatible with the expected values, be- 
cause the EW range from 71 eV (in the QSO subsample) to 283 
eV (in the z<l subsample). 
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Fig. 10. Fit with the 'two-component model'. Top panel: sub- Fig. 11. Fit with the 'two-component model'. Top panel: Sey 
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In addition, the EW estimated with our simulations are gen- 
erally higher when using a broader spectral window, which is 
compatible with the presence of a broad emission line. 



5.2. Discussion based on phenomenoiogicai models 

We considered a simple phenomenoiogicai Gaussian to fit the 
iron line and calculated the EW considering a monochromatic 
continuum at 6.4 keV. In these fits, we found that it is required 
by the data with high significance: from the spectral fitting to a 
simple Gaussian, the improvement of adding the Gaussian to the 
continuum is at 1 - P(Av,Ax^) - 9 x 10"'^ in the full sample 
(see Table|3]l, value that corresponds to a 6.8cr confidence level. 
The EW for the narrow line, is 95+22 eV for the full sample and 
lays between 77+35 eV (for the Abs subsample) and 167+67 eV 
(for the z< 1 subsample): again, such values are consistent with 
those of fluorescent lines from material far away from the central 
BH. 

Allowing for a broad neutral line {Efixed-crfree second row of 
each subsample in Table[3]l, the errors on the line physical width 
could not be constrained in the full sample and in the majority 
of its subsamples: the upper limit of ~ 300 eV found in the full 
sample indicates that we cannot exclude the presence a broad 
line component in our detected lines although we are not able to 
constrain its properties. Looking at Table |3] Col. 4, the subsam- 
ples where a significance higher than 3cr is found in the fit with 
E fixed - o'free (sccond row of cach subsample) are the Abs sub- 



sample and QSO subsample. In both cases, there is a degeneracy 
between the iron line and the reflection component. 

The EW values obtained from the fits with Efj^-o-free gener- 
ally are in agreement with those found in prev ious works based 
on averaging methods, e.g. lCorral et al.l(l2008h . being of 107+50 
eV in the full sample. A minimum of 36 % of this value is in 
principle the contribution from lines other than the K„ neutral 
iron line, such as the Fe K/3, the Ni K,,, the Compton Shoulder 
of the Fe Kq, (Matt 2002), and that is expected to contribute to 
the detected line. The EW of 107+50 eV estimated for the full 
sample would mean that the contribution of the K„ iron line only 
is EW ~69 eV. We expect EW encompassing a range of values 
from about 50 eV to 100 eV for non-obscured AGN, including 
the measurements presented here. 

Investigating subsamples in luminosity and redshift, we 
found a hint of a difference in the line EW between the Seyferts 
subsample and the QSOs one: the line EW, for example in the 
second line of each sample in Table[3] is 88+32 eV in QSOs and 
292+102 in Seyferts. Thus, the EW decreases with the increas- 
ing luminosity of the X-ray continuum in agreement with previ- 
ous results f or narrow ( Iwasawa & Taniguchi 1993) and broad 
lines (iJimenez-Bailon et al.i 2005). A similar result is found also 
from the model-independent estimates of the EW presented in 
the third column of Table |2] 

In most cases, if the line width is left free to vary in the spec- 
tral fitting, its best fit value becomes higher than zero, as dis- 
cussed above; on the contrary, if the centroid energy is left free, 
its output value does not differ from the input 6.4 keV energy 



10 



S. Falocco et al.: Iron lines in tlie XMM CDFS deep survey 



(see second and third lines of Table [3]l. This indicates that the 
majority of the line flux is emitted for fluorescence from neutral 
iron. 

The EW values (e.g. focusing on the Eji^ed — o- fixed fitting 
presented in Col 9 of Table |3]l are consistent with the ones cal- 
culated making use of the simulations, presented in Table |2] and 
explained in Sect. 14. II this confirms that our estimates are robust. 

The EW measured in this work are broadly consistent with 
the values of individual bright sources in the literature: this con- 
firms that the iron lines are emitted by the majority of the sources 
in our sample because the EW of the minority Iron-line-emitting 
sources would be much higher than that observed in individual 
sources. Our resulting EW are also consi stent with previous es- 
timates made in iron Hne stacking, e.g. in lCorral et al.l(l2008h . 

5.3. Discussion based on piiysical models 

We used a model more specific than a simple Gaussian to fit 
our data and calculated the EW considering a monochromatic 
continuum. The 'physical' model is composed by a narrow core 
centred at 6.4 keV and a relativistic component (diskline in 
xspec). 

This two-component model is expected in AGN, where re- 
flection of the primary emission is expected to occur both in axis- 
symmetric regions far away from the BH (such as the torus) and 
close enough to it for the detection of an y line broadening, e.g. 
in the accretion disk dFabian et al.lll989h . Both components are 
in principle expected. Our values of the EW reported in Table 
|4] support this, being consistent with previous results based on 
averaging techniques and on the study of individual higher qual- 
ity spectra of AGN in the local Universe. Moreover, our model- 
independent estimates between 173 and 236 eV (full sample, 
bulk range in Table |2]i are inconsistent with having only a nar- 
row component, which does not commonly exceed 100 eV in the 
local Universe. This provides an interesting indication, indepen- 
dent of the spectral fitting, of the presence of a relativistic line 
component in the XMM CDFS data. However, the improvement 
of the fit by including the relativistic component (Col. 3 in Table 
|4|i is not statistically significant (e.g., for the full sample, 2.1 cr 
confidence level), so we have not formally detected a broad com- 
ponent our samples. Having said that, the two line components 
are strongly correlated and we cannot really disentangle them: 
we can see the 1, 2, 3 cr contours (for two interesting parame- 
ters) for the full sample in the space formed by the EW of the 
relativistic line and that of the narrow line in Fig. [T2](top panel). 
There are equally-good fits including only each of the compo- 
nents or a combination of both for most samples (see cols. 1, 2 
and9ofTable|4]i. 

It is important to notice that our contours (shown for the full 
sample) exclude the point with zero EW simultaneously for both 
components at much more than three cr confidence level. This 
indicates a very secure detection of an iron line feature. 

For the full sample, the EW of the narrow line with no rela- 
tivistic component is 95 + 22 eV (see first row of Table |3]l while 
the EW of the relativistic line fixing the flux of the narrow line to 
zero is 347 + 105 eV (Col. 10 in Tableg] first row). These results 
are, within the uncertainties , compa tible with the measurements 
and uppe r limits found in ICorral et al. (2008) and Chaudhary 
et al. (12012) . Taking into account that our values include con- 
tributions of fluorescence in material close to the central SMBH 
and in outer regions that we are not able to disentangle, they are 
also consistent, considering the error bars, with the "average" 
EW=76-1 43 eV of the relativist i c fine s detected in individual 
sources in lde la Calle Perez et al.l(l2010l) . 




50 100 150 

EW of narrow component (eV) 



200 




Fig. 12. Contours at 1 (continuous line), 2 (dashed line), 3 (dot- 
ted line) sigma of the EW of the narrow and the relativistic line 
in the 'two-component model' (top panel). Same contours for 
the reflection factor R of both the non-relativistic and the rela- 
tivistic components of the comprehensive model (bottom panel). 



5.4. Discussion based on compreiiensive models 

We fitted the average spectra with a more complete model with 
two reflection components. The first one is emitted from the in- 
nermost regions of AGN and smeared by relativistic effects due 
to the vicinity of the central SMBH, the second component is 
instead emitted from outer regions and is un-blurred. The re- 
flection factor contours at 1, 2, 3 sigma, are shown in Fig. [12] 
bottom panel. The figure shows that the relativistic component 
is different from zero at two sigma confidence level, and the 
non-relativistic one at three sigma. The two components are cou- 
pled: like the above described 'two component' model of the 
iron line, these cannot be accurately disentangled. Anyway, we 
should note that the point with zero reflection factor for both 
components is excluded at much more than three sigma confi- 
dence level, indicating a secure, highly significant detection of 
reflection features. 
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6. Conclusions 

We studied a representative sample of 54 AGN (100 spectra from 
XMM-Newton) with the best spectral S/N in the ultra-deep ob- 
servation of the Chandra Deep Field South. We investigated sev- 
eral subsamples: two subsamples of column density of the intrin- 
sic absorber (borderline at log (A^h) - 21.5), two subsamples in 
luminosity (with cut at log (L) = 44), and finally two in redshift 
(segregated at z = 1). 

We averaged the spectra of the full sample and subsamples 
deve lopin g further our own rnethod previously used in Corral 
et al. (l2008i) and iFalocco et al.l (l2012h to characterise the instru- 
mental and averaging effects. To do this we made use of com- 
prehensive simulations. 

With our new methodology and best data quality so far we 
could find convincing and solid evidence, using X-ray stacking, 
for a highly significant iron line in distant AGN, taking into full 
account the instrumental and methodological effects. In the spec- 
tral fits we convolved our models with a gaussian smoothing 
with a variable !,{£) introduced by the limited spectral resolu- 
tion of the X-ray instruments and the averaging process. 

We can summarise our conclusions as follows: 

- The iron feature is significantly detected in the full sample 
and in all its subsamples. 

- The iron line, modelled with an unresolved Gaussian at 
6.4 keV (as expected from neutral material far from the 
central BH), is detected at 6.8cr in the full sample, with 
EW=95 ± 22 eV. In the subsamples, the EW measurement 
ranges from 77 eV for Abs (our objects with log (Nh) > 
21.5) and 167 eV for the sources with z < 1 respectively. We 
find evidence for the Iwasawa-Taniguchi effect comparing 
the EW of the low luminosity (Seyfert) and high luminos- 
ity (QSO) samples. Allowing for a resolved Gaussian profile 
improves the fit slightly, but not significantly ($ 2o"). 

- Using a combination of an unresolved line from neutral iron 
as above and a relativistic line, to take reflection in mate- 
rial far and close to the central BH into account, we found 
again that the broad component, despite improving the fit 
moderately in most cases, is not formally required. However, 
we also found that the fitted parameters of both components 
are strongly coupled, finding good fits with either or both 
components. We cannot really disentangle both contributions 
with our current data. Using a single relativistic component 
with no unresolved line, we find EW= 347 + 105 eV for the 
full sample. 

- The E W of the narrow iron line and that of our tentative rela- 
tivistic line are both broadly consistent with those (including 
the upper limits) estimated in previous works ba sed on aver- 
aging and stacking methods (Corral e t al.ll2008 : Chaudharv 
et al. 2012) and in individual bright sources in the literature 
(Ide la Calle Perez et al.ll2010 ). which leads to the conclu- 
sion that an iron feature must be widespread among spectra 
of AGN. 

- The model-independent EW estimates made using our thor- 
ough simulations produce similar values as those from model 
fitting described above, confirming the robustness of the re- 
sults presented in this paper Interestingly, the EW is incon- 
sistent with having a narrow line only and becomes higher 
when a broader energy range around the line centroid is con- 
sidered. This suggests the presence of a broad iron feature. 

Concluding, although the observed lines are well modelled 
as fluorescent lines at 6.4 keV in material far away from the 
primary emission, we found also some tentative evidence for a 



component emitted relatively close to the SMBH. Our work has 
confirmed that high S/N and a large amount of spectral counts 
are required to characterise average iron features in distant AGN. 
For a more detailed characterisation of its physical origin, it is 
necessary to accumulate a higher number of spectral counts in 
the average spectrum. 

Although the complexity of the X-ray AGN spectra often 
makes challenging the determination of the properties of iron 
features, with the deepest XMM-Newton data available and a 
careful characterisation of methodological effects, we have con- 
firmed that iron emission lines represent a universal feature of 
AGN in the distant Universe. 
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Table 3. Results of fits of the average spectrum of the full sample and its subsamples with gsmooth(8i(phaxpex+gauss). 



Sample 


x-ly(c) 


x^ly(g) 


l-P{hv,hx-) 


Nh 
10--cm" 


(1) 


(2) 


(3) 


(4) 


(5) 



9 X 10"'^ (6.80-) 
[0.4 (0.9o-)] 
[1 (Oo-)] 
[0.6 (0.52O-)] 



36+"" 

0.36 + 0.22 
0.36 + 0.21 
0.36 + 0.21 



(6) 

1 63+"'*-' 

1 63+"" 
1.63 + 0.06 
1.63 + 0.05 



SJ) 

1 00+" ™ 
1.0+3* 
1 0+''-' 

1 0+3 ' 

'■ -in 



E 

keV 

(8) 



eV 
(9) 



EW 


<<-> 


<log(L)) 


eV 




erg s"' 


(10) 


(11) 


(12) 



Full 133.4/95 



87.0/94 
86.2/93 
87.0/93 
85.9/92 



=6.4 


=0 


=6.4 


09+236 


6.4+0.20 


= 



3 X 10-" (5.50-) 
[1 (Oo-)] 
[1 (Oo-)] 
[1 (Oo-)] 



6 44+"°* 



91 



-197 

-«7 



95+22 
107+50 
94+22 
109 + 44 



1.34 43.74 



Unabs 139.2/95 



108.7/94 
108.7/93 
108.7/93 
108.7/92 



< 0.82 
<0.82 
<0.82 

< 0.82 



1.87 + 0.03 
1.87 + 0.03 
1.87 + 0.03 
1.88+0.03 



1.0 + 0.4 



1 0+"-^ 

'■"-'-0.2 
1+0.2 



=6.4 
=6.4 

6 4+"-3 



=0 

1+822 
'-1 

=0 

^+800 



107 + 32 
107 + 32 
107 + 32 
107 + 32 



1.27 43.78 



Abs 



109.4/95 



94.2/94 
85.4/93 
94.2/93 
85.6/92 



1 X lO"" (3.90-) 
[3 X 10-3 (3o-)] 
[1 (Oo-)] 
[3 X 10-3 (2.9o-)] 



3.67+0.55 
3.39+0.53 
3.67 + 0.55 
3.42 + 0.54 



1.49 

1.39: 



t0.08 

:0.08 



I co+O.IO 
-^ -0.09 

1.39 + 0.08 



1+..6 

0.003!i;™ 

1+3.7 

0.003^^ 



=6.4 
=6.4 
6.4 + 0.2 
6.47 + 0.20 



=0 

455+-"'' 
=0 



77 + 34 1.28 43.74 

243 + 123 
77 + 32 
239 + 109 



z<l 



112.6/95 



97.0/94 
94.0/93 
94.0/93 
94.0/92 



8 X 10-' (3.90-) 
[0.08(1.70-)] 
[0.08(1.70-)] 
[0.22 (1.2o-)] 



0.86 + 0.40 
0.92 + 0.40 
0.92 + 0.40 
0.92 + 0.40 



1.58 + 0.08 
1.56 + 0.08 
1.56 + 0.08 
1.56 + 0.08 



1 0+'3-5 
'■"-1.0 

'■"-0.6 

'■"-0.6 

1.0+S' 



=6.4 
=6.4 

6^4o±;^-?^ 



=0 

4+150 

=0 

2+?-- 



167+67 
81+34 
81+34 
81 + 34 



0.64 43.21 



z>l 



109.9/95 



78.5/94 
76.2/94 
78.5/93 
76.2/92 



2 X 10-" (5.60-) 
[0.1 (1.5cr)] 
[1 (Oo-)] 
[0.32 (Itr)] 



< 0.24 

< 0.32 

< 0.441 

< 0.32 



1.63: 
1.60: 



0.06 
0.05 



2 X 10-'-' (60-) 
[0.066(1.840-)] 
[1 (Oo-)] 
[0.17(1.40-)] 



36-"-3" 

"■^"-0.29 

<0.47 
36+"-' 

"■^"-0.29 

<0.56 



1 63+"" 

'■"^-0.04 

1 60+""' 

'■""-0 04 



0:008+i;™ 
3+17 

0.008!ii 



TT^ 

'■"-0.8 

004+'""° 

"■"" -0.004 

1 0+^ = 

'■"-0.8 

004+'""" 

"■""^ -0 004 



=6.4 
=6.4 

g40+0.25 

6.40 



,:8:?5 

07 



=0 



=0 

151+"^ 
'■^'-iM 



114+30 
114+42 
150+40 
150 + 55 



1.80 44.07 



Sey 



129.9/95 



93.9/94 
90.6/93 
93.9/93 
90.4/92 



1.58 + 0.07 
1.51 +0.07 
1.58 + 0.07 
1.51 +0.06 



Twrr 



1 60+""' 

'■""-0.02 

1.65 + 0.04 
1.68 + 0.4 
1.65 + 0.04 



14+1-* 
"■'^-0.1. 

0.66+;-3' 

0.9+'3 

0.61 



=6.4 
=6.4 
6 40+"°' 

"■^"-0.05 

6 5+"°' 

" "-^-0 70 



=0 

450+ 

=0 

425J 



187 



103+27 0.91 43.35 

293+102 

103+27 

289+96 



QSO 204.2/95 



178.2/94 
166.6/93 
167.3/93 
164.9/92 



3x 10-' (5.IO0-) 
[5.3 X 10-* (3.5o-)] 
[9 X 10-* (3.3o-)] 
[0.0013 (3.220-)] 



0.11 



^-0.08 

0.20 + 0.15 
0.24 + 0.15 
0.19 + 0.15 



-0.04 



=6.4 
=6.4 
6.40+"-°5 



-0.30 



6.34 



;8:i? 



=0 

ioT+120 
'"'-100 
=0 

130+?' 



119+30 

88+32 
75+22 
93 + 31 



2.02 44.34 



Notes. Columns: (1) sample; (2) ^^/dof of the fit with the continuum; (3) ^^/dof of the fit with the continuum and the Gaussian; (4) probability 
P{iS.x^, Ay) (of adding the Gaussian to the model and of leaving each parameter free to vary, see text); (5) intrinsic column density estimated in the 
fits; (6) slope of the powerlaw; (7) reflection factor; (8) central energy of the Gaussian; (9) a of the Gaussian; (10) EW of the Gaussian; columns 
1 1 (average redshift of the sample) and 12 (average logarithmic luminosity of the sample in erg s-'), are repeated from Table 1 for reference 

Table 4. Fits results of the average spectrum of the full sample and subsamples using gsmooth®[phax(pex+laor)+gaus)]. 



Sample 
(0) 


(1) 


(2) 


i-P(Av,Ax^) 
(3) 


Nh 

lO^^cm-^ 

(4) 


F 

(5) 


Rel„fi 
(6) 


EW„, 

eV 

(7) 


EW„„ 

eV 

(8) 


2 / 

y Jv 

■^une-rc'l' 

(9) 


EW„„„_„, 

eV 

(10) 


Full 


87.0/94 


82.6/93 


3.6xl0--(2.1o-) 


0.24 ± 0.24 


1.61 ±0.06 


1 0+*'-' 


140 ± 120 


67 ±28 


98.4/94 


347 ± 105 


Unabs 


108.7/94 


108.7/93 


l(Ocr) 


<0.82 


1.87 ±0.03 


1+0.4 


n+ioou 


107:- 


138.9/94 


n+io(K) 


Abs 


94.2/94 


92.6/93 


0.2(1.3o-) 


2.84 ± 0.60 


1.29 ±0.08 


0-^i 


404 ± 140 


< 1000 


92.6/94 


404 ± 140 


Z<1 


97.0/94 


95.2/93 


0.2(1.3o-) 


0.68 ± 0.40 


1.55 ±0.09 


i:: 


173 ± 173 


104 ± 93 


99.7/94 


362 ±210 


Z>1 


78.5/94 


n.6/93 


0.3(lo-) 


<0.3 


1.58 ±0.04 


0.007!--' 


84 ±84 


98-; 


94.6/94 


318± 119 


Sey 


93.9/94 


92.2/93 


0.2(1. 3o-) 


o-26:i|-^^ 


1.57 ±0.07 


i:? 


118± 118 


82-= 


107.6/94 


380 ± 130 


QSO 


178.2/94 


166.7/93 


7xl0-''(3.4o-) 


< 0.138 


1.58 ±0.08 


0.13t 


228 ± 70 


49 ±45 


169.7/94 


316 ±82 



Notes. In all fits we have fixed R; =6,/? =-2, 9 =30. Columns: (0) sample; (l);^'^/dof of the fit with the Gaussian with Efi„d - (rfi„d only, repeated 
from Table|3]for reference; (2) ^^/dof with the two-component model (see text); (3) Significance of the laor (expressed in l-P(Av, Ax^) and in 
(t); (4) intrinsic column density of the neutral absorber. (5) powerlaw slope. (6) reflection factor; (7) laor EW measured from the two-component 
model; (8) EW of the narrow component; (9) ^'/dof of the fit after having fixed the flux of the narrow line component to zero; (10) laor EW 
measured from the same fit of Col. 9. 
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Table 5. Fits results of the average spectrum of the full sample and subsamples using gsmooth®[ pha x kdblur2 
®(pexm)+pexm)] . 



Sample 


x^ly 


Nh 


F 


^non-rel 


Rrd 


(0) 


(1) 


(2) 


(3) 


(4) 


(5) 


Full 


84.1/94 


0.46±0.03 


1.71+0.06 


-0.35+0.18 


45+*"** 


Unabs 


104.0/94 


<0.10 


1.95+0.06 


-0.50+0.27 


43+0.54 


Abs 


90.2/94 


3.68 ±0.53 


1.56+0.09 


-0 167+"'*** 

" -0 228 


59+0:51 

0.08!i|-* 


z< 1 


91.5/94 


1.01+0.40 


1.65+0.10 


-0.317+0.220 


z> 1 


80.5/94 


0.343±0.290 


1.76+0.07 


-0.53+0.23 


Sey 


89.0/94 


0.48+0.30 


1.69+0.09 


-0.39+0.22 


0.54 + 0.45 


QSO 


73.9/94 


0.57+0.32 


1.76+0.08 


-0.33+0.26 


1 6+"'*' 



Notes. In the fits we have fixed: Ri„ = l.2'i5Rg, Rbr=20Rg, Roiu=400Rg. Columns: (0): sample; (1): x /^ of the fit; (2): intrinsic column density in 
lO^^cra"^; (3): photon index; (4): reflection factor of the reflection component produced far away from the central engine; (5): reflection factor of 
the reflection component produced close to the central engine. 
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